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I. INTRODUCTION 



Landau levels are equally spaced, highly degenerate energy levels of a charged parti- 
cle moving on a plane under the influence of a perpendicular uniform and static magnètic 
field [1,2]. They play very important roles in our current understanding of two dimensional 
electron systems and their characteristic behaviors such as quantum Hall effect [2-4] and 
superconductivity [5] and in exploring electrònic properties of various bound states of mat- 
ter such as àtoms, molècules and condensed matter in strong magnètic fields [6]. Regarding 
recent experiments [7,8] devised to measure Wigner function of various states of light and 
matter and to observe its remarkable properties, determination of Wigner functions of Lan- 
dau levels and search for their salient features become very promising since they can be 
realized in several classes of systems. 

Wigner function is quantum mechanical analogue of classical phase-space probability dis- 
tribution function. It is the central concept of Weyl-Wigner-Groenewold-Moyal (WWGM) 
quantization which was initiated just after the formulation of quantum mechanics [9]. 
WWGM quantization consists of well established correspondence rules between the opera- 
tor formulation of quantum mechanics and its phase-space formulation [9-11]. According to 
these rules Wigner function can be computed from 

W X y = J n Jx^ + \y)^A^-\y)^dV{y), (1.1) 

where, H being the Hamiltonian operator, %p\ is an eigenfunction of the Schròdinger equation 
Hip\ = E\ip\, ip denotes the complex conjugation of ip and (q, p) = (çi, çd,Pi, ■■■,Pd) 
are the phase space coordinates. dV(y) stands for D-dimensional volume element in in- 
variables, h = h/2n is the Planck's constant, and y ■ p denotes the usual scalar product of 
vectors y and p. The functions W\y = W\\i(q, p) are called off-diagonal Wigner functions 
and what is referred to as Wigner function is known also as the diagonal (or, pure state) 
Wigner function which corresponds taking ip\ = ifj\> in (1.1), where À may be a multi index 
set of quantum numbers. Throughout this paper the pair (qj,Pj) will denote the canonically 
conjugate position and momentum variables for each j = 1,2,..., D, bold face letters will 
stand for vector quantities and when there is no risk of confusion arguments of functions 
will be suppressed. 

Deformation quantization which was formulated in the pioneering work [12] in the sev- 
enties by generalizing WWGM quantization has become the third autonomous and logically 
complete formulation of quantum mechanics beyond the conventional ones based on opera- 
tors in Hilbert space or path integral [13,14]. In this formulation quantum effects are encoded 
in a new operation called *-product. This is a noncommutative but associative composi- 
tion rule between functions defined on phase-space which can be any symplectic manifold 
[14,15] (see also [13,16] for recent developments) . In terms of *-product the spectrum and 
corresponding phase-space eigenfunctions of a Hamiltonian function H can be determined 
through two-sided *-eigenvalue equation [17-19] 

H*W X = W X *H = E X W X . (1.2) 

In that case the Wigner function of H corresponding to the eigenvalue E x is the *- 
eigenfunction W\ of (1.1). When the underlying phase-space is R 2D with globally defined 
coordinates (q, p), the *-product is given by 
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* = exp[-^^(9 g ,9 P , - (1.3) 

where we use the abbreviation d Xi = d/dxi and the convention that 9 and d are acting, 
respectively, on the left and on the right. (1.3) is known as the Moyal *-product. 

Apart from a few simple cases the use of (1.1) to obtain Wigner function is rather difficult 
since for a given wavefunction the resulting integrals are not easy to cope with. On the other 
hand, only a few *-eigenvalue problems have been analytically handled until now. In fact, 
the bàsic physical systems such as H atom and harmònic oscillator were considered in [12] 
with special emphasis on determination of their spectra. Later on, some studies utilizing 
different techniques considered the same and some other systems in phase-space formulation 
[17-21]. To the best of our knowledge, the Ref. [21] is the first in dealing with Wigner 
function of Landau levels by using some WWGM correspondence rules. Systematic use of 
*-eigenvalue equations were taken up in Ref. [19] which inspired us much. 

The main goal of the present paper is to generate all Wigner functions and corresponding 
marginal probability densities of Landau levels on all possible phase-space planes. In doing 
that we shall follow an autonomous approach by working entirely within the framework of 
deformation quantization. All genèric properties such as reality, normalization, projection 
and orthogonality as well as transformation and symmetry properties of these bound state 
Wigner functions are rigorously investigated in this framework without any reference to 
wavefunctions or to WWGM correspondence rules. After solving a *-eigenvalue equation for 
the ground state Wigner function, Wigner functions of the excited states will be produced 
algebraically by successive applications of creation and annihilation functions. Then, off- 
diagonal Wigner functions and marginal probability densities will be obtained by means of 
associated generating functions. 

The organization of the paper is as follows. In the next section the main points of defor- 
mation quantization, mainly those needed in subsequent investigation are briefly reviewed. 
Formulation of the problem and a closed algebraic form of diagonal Wigner functions of 
Landau levels are given in Sec. III, and their explicit functional forms are derived in Secs. 
IV and V. Symmetries and other genèric properties are taken up in Secs. VI and VII. A 
generating function is introduced in Sec. VIII and all (diagonal and off-diagonal) Wigner 
functions are generated by means of it. In Sec. IX phase-space coherent state property of 
this generating function is presented and it is shown that its integrated forms serve as gener- 
ating functions for marginal probability densities derived in Sec. X and presented altogether 
in Tables I and II. Final section is devoted to some remarks concerning justifications of some 
our results, implementation of unitary similarity transformations and gauge transformations 
in a phase space and to some general implications for Wigner function theory. In appendix 
we derive an operator identity used in the main text and exhibit a new finite sum formula 
for the generalized Laguerre polynomials. 

II. STAR PRODUCT AND MOYAL BRACKET 

Let us denote the linear space of complex valued, smooth (differentiable to all orders) 
functions defined on R 2D by Af. The Moyal *-product (1.3) is a bilinear product * : AfxAf — > 
Af, which is, for all fj G Af, associative 
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h * (fi * fs) = (h * h) * fs, 



(2.1) 



and obeys the so-called closedness property [17,22] 

L * * í* dv = L f^ dv = L h * hdv - (2 - 2) 

This can be proved via integration by parts and it implies that integration, all over the 
phase-space, is a trace on the *-algebra of functions. From (1.3) immediately follows that 
under complex conjugation the *-product of functions satisfies the following relation 



(/i*/2) = / 2 */i- (2-3) 

It is also evident that /i(q) */ 2 (q) = /i(q)/2(q) and ^(p) *# 2 (p) = 0i(p)02(p), where the 
product on the right hand sides is the usual commutative pointwise product of functions. 
In terms of *-product Moyal bracket {, }m is defined, for all /, g e jV, as 

{f,d}M = f*9-9*f. (2.4) 

The properties of the ^-product ensure that Moyal bracket (MB) is bilinear, antisymmetric 
and obey the Jacobi identity and Leibniz rule. With respect to MB, Af acquires a Lie 
àlgebra structure. The most important properties of ★-product and MB are the following 
limit relations 

]imf*g = fg, (2.5a) 

fi — >0 

hmUh)- 1 {f,g} M = {f,g} P , (2.5b) 

n— >U 

where 

D 

{f,g}p = J2( d iJ d p k 9 - d Pk fd qk g), (2.6) 

k=l 

is the usual Poisson bracket (PB) of the classical mechanics. The limit relations (2.5) reveal 
the fact that the associative ^-àlgebra and Lie àlgebra structure of H given by MB are, 
respectively, deformations (in the sense of Gerstanhaber [23]) of associative àlgebra struc- 
ture of Aí with respect to above mentioned pointwise product and of Lie àlgebra structure 
determined with respect to PB. 

The set A of phase-space functions /'s which satisfy 

{f,g} M = ih{f,g} P , (2.7) 

for all g e H plays an important role in deformation quantization. First, A is a Lie 
subalgebra with respect to PB which is preserved under the deformation PB — > MB. 
Second, with respect to PB, A acts as a derivation on the ^-àlgebra of functions, that is, 

{/, 9i * 92} p = {/, 9i}p *92 + 9i* {/, 92} p, (2.8) 

for all g±,g2 & Af and / e A. Finally, quantum mechanical time evolution of any element of 
A is a classical evolution. This implies that when the Hamiltonian function belongs to A, 
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the classical and quantum time evolutions of an observable coincide. In the nomenclature 
of the deformation quantization A is called the invariance àlgebra of the *-product and its 
elements are called preferred (good or distinguished) observables. As can easily be verified 
from (1.3) and (2.4) A is spanned by 

{1, Qj, Pj, QjQk>j, PjPk>j, QjPk}- (2.9) 

These define 2D 2 + 3D + 1 dimensional affine symplectic àlgebra wd © sp(2D), where wd 
stands for 2D + 1 dimensional Heisenberg-Weyl àlgebra and sp(2D) denotes 2D 2 + D di- 
mensional symplectic àlgebra. 

For g — Ci • p + /, where Ci is a constant vector and / = /(q), we have 

00 1 ifi D «- -» 
9*9 n = 9 n+1 + E 7f( 9")V(E d qk d Pk ) j 9 n . 

j=2 J- Z k=l 

Hence, if / is first order in q^s then g * g n = g n+1 for any n, which implies that 

( 9ir ) n = g n , for 3 = cq + C! -P + C2 -q, (2.10) 

where co is a constant, c 2 is another constant vector and we have defined the *-power of 
g e M as 

(^*) n = ^*^*···*^, (ntimes). (2.11) 

Finally in this section we should note that for two functions of the form / = ai + h±, g — 
a 2 + h 2 such that a±, a 2 G A, and {/ii, /í2}m = ih{hi, h 2 }p we have {/, g}u = ih{f, g}p- As 
particular cases this holds for hj = hj(p), or for hj = hj(q). 



III. WIGNER FUNCTIONS OF LANDAU LEVELS 

The well known Landau Hamiltonian H L for a spinless particle of charge q > and mass 
m moving on the gig 2 -plane is (in the Gaussian units) 

Hl = ^(P-f A ) 2 = ^(«i 2 + «2 2 ), (3-1) 
where c is the speed of light, A = A(q) is the vector potential of the magnètic field 

B = d qi A 2 - d q2 A u (3.2) 
and v = (p — |A)/m is the velocity vector whose components obey 

{vi,v 2 } M = i—^-B. (3.3) 

In writing (3.3) we have used the fact that the phase-space of the Landau problem is R 4 
with canonical coordinates (q, p). In that case *-product and MB are given by (1.3) and 
(2.4) for D = 2. We should also note that (3.3) is vàlid for any B. 
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When the magnètic field is constant the general solution of (3.2) is 

A = (~Bq 2 + d qiX , \B qi + d q2 x), (3.4) 

where x = x(q) i s an arbitrary gauge function. In such a case we have, in any gauge \i t wo 
constants of motion 

Xi = m(v 2 + ooqi), (3.5a) 
X 2 = -m{vi - ujq 2 ), (3.5b) 

where u = qB/mc is the cyclotron frequency. (Xi/mcu, X 2 /mcj) correspond to the coordi- 
nates of cyclotron centre and satisfy the following gauge-independent relations 

{Xi, X 2 }m = — imhuj, (3.6a) 
{ Vj , X k } M = = {H L , X k } M , j, k = 1,2. (3.6b) 

The relations {H L , X k } M = imply that H L ,X 1 and X 2 are constants of motion and since 

dXi A dX 2 A dH L = m 3 uj[(v 1 dq 2 — v 2 dqi) A dv x A dv 2 + — v • dp A dq\ A dq 2 ] (3.7) 

they are functionally independent provided that v ^ 0. In Eq. (3.7) dqj,dpj denote the 
coordinate differentials (coordinate 1-forms) and A stands for the antisymmetric exterior 
product [15]. The above observations emphasize the quantum superintegrability of the 
problem. Moreover, by the remark at the end of previous section, it is easy to check that 
MBs in Eqs. (3.3) and (3.6) are equal to %h times the corresponding Poisson brackets. This 
implies that {Hl, X±, X 2 } are at the same time classical constants of motion in any gauge, 
and hence the problem is also classically superintegrable. As a result, the Landau problem 
is one of the rare problems for which the classical and quantum superintegrability coincide. 
For the phase-space quantization of the Landau problem we shall use Hl and 

J = —(X 2 -m 2 v 2 ), (3.8) 

as a complete Moyal-commuting phase-space functions. J corresponds to canonical angular 
momentum q\p 2 —q 2 pi in the symmetric gauge x =constant. We then introduce two mutually 
commuting pairs of dimensionless creation and annihilation functions 

1 / 

a = — (ui + tv 2 ), 

à = — (i>i — iv 2 ), (3.9a) 
70; 

b = —(X 2 + tX 1 ), 

b = —(X2-ÍX!), (3.9b) 

where 7 = (2hc/ 'qB) l l 2 = (2h/muj) 1 ^ 2 is the so-called magnètic length. 7 corresponds the 
radius of a disc from which a flux quantum hc/e passes for q = e, e being the elementary 
charge. The nonvanishing Moyal brackets of (3.9) are 
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{a,a} M = 1 = {b,b} M . (3.10) 
In terms of real number functions defined by 

N a = a*a, N b = b*b, (3.11) 

(3.1) and (3.8) take the form 

H L = hu(N a + ^), (3.12a) 

J = h(N b -N a ). (3.12b) 



By making use of 



{a,(a*) k } M = k(à it ) k - 1 , (3.13a) 
{b,(h) k } M = k(K) k -\ (3.13b) 



and by defining the ground state Wigner function W by 

a*W = = b* W , (3.14) 
it is straightforward to check that 

W nl = 4ï(<i*) n * (6*)' * * (a*) n * (6*)', (3.15) 

77/. £ . 

satisfy the following "*-ladder structure" 

a*W nl = Wn-^i + a, 

6*W ní = W niI _i*6, (3.16a) 
a * W n ; = W^+i^ * a, 

b*W nl = W nJl+1 *b, (3.16b) 

where n, / are positive integers. We assume that Wo = Woo is a real function. Eqs. (3.16b) 
can also be verified by taking the complex conjugates of (3.16a) and then performing the 
shifts n — > n + 1 and l — > Z + 1. For the number functions (3.11) we have 

iV a * W„i = W n j * N a = a* W n -!,i *a = nW nh (3.17a) 
iV fe * W nl = W nl *N b = b* W^-i * 6 = ZW ní . (3.17b) 

Hence, {W n \ : n, l — 0, 1, ...} is the set of simultaneous eigenfunctions of iV a and iV 6 whose 
spectra are bounded from below. 

In view of the above relations it is now easy to show that W n i satisfy 

H L * W nl = W nl *H L = E n W nh (3.18a) 

J * W nl = W ni * J = JnlW nh (3.18b) 

with *-eigenvalues 
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E n = hu(n+^), (3.19a) 
J nl = h(l-n), n,/ = 0,1,2, ... (3.19b) 

E^s are the well-known infinitely degenerate ( since every E n is independent from the 
quantum number /) Landau levels. The corresponding Wigner functions are given by (3.15) 
(or, by Eq. (3.22) below) in a closed algebraic way and in a gauge-independent manner. We 
should note that the degeneracy of a Landau level is a reflection of the fact that energy is 
independent of the location of the guiding center of gyrating charge. It is infmite when the 
motion takes place on the whole R 2 plane and otherwise it is finite. 

In the next section we will choose the symmetric gauge x =constant, that is 

(A 1 ,A 2 ) = ^(-q 2 ,q 1 ), (3.20) 

and discuss the Wigner functions in any gauge in section XI. Before doing that we should 
note that in a gauge which is at most quadratic function of q/s a,a,b,b are some linear 
functions of q/s and p/s. Hence, by virtue of (2.10) we get 

(à*) k = à k , (a*) fc = a k , (3.21a) 
(&*)* = b k , (bA k = b k , (3.21b) 



and rewrite (3.15) as 



W ni = — -a n *b l *W *a n *b l . (3.22) 
n\l\ 



Note that the *-product between a's and ò's can be removed for a * b = ab and a-kb = ab. 

IV. THE GROUND STATE WIGNER FUNCTION 

Let us introduce the complex coordinates 

z = -^(çi +iq 2 ), P= ^(Pi-m), (4.1a) 

z = -^={q 1 -iq 2 ), p= -^={p! + ip 2 ), (4.1b) 



which imply 



d qi = -^( d * + d z), d Pl = —(d p + d p ), (4.2a) 







+ d P ) 


d P2 = 


7=2 {dp 





d q2 = ^(d z -d s ), d P2 = -j=(d p -d p ). (4.2b) 
In terms of new coordinates the *-product (1.3) (for D = 2) transforms to 

* = exp[Uh(dzd P + dzdp - d P d z - dpdz)}, (4.3) 
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and the annihilation functions take the form 



\[2 imu 

:(p ( 4 - 4a ) 



mju 



In that case the defining relations (3.14) of W are as follow 



\[2 , imu . /, \ 

-— z). (4.4b) 



p*W o -—z*W o = 0, (4.5a) 

2 777 Cl? 

p*W -— -z*W o = 0. (4.5b) 
Combining the complex conjugate of (4.5a) with (4.5b) we get 

711 (jj A 

d z W = -—zW , d p W = -pW . (4.6) 

n mun 

By ignoring a possible factor f(z,p) which destroys the reality of Wq, and by defining 

„ 1 1 2 p 2 mu 2 2 , À _. 
H = -PP + -mu 2 zz = + ^q 2 , (4.7) 
m 4 2m 8 

which corresponds the Hamiltonian function of a two-dimensional isotropic harmònic oscil- 
lator, the general real solution of Eqs. (4.6) is 

W = 4e'^ H °. (4.8) 

We have normalized Wq as J R4 W dV = ti 2 , where the volume form dV is 

dV = dqi A dq 2 A dp 1 A rfp 2 — dz Adz Adp A dp. (4.9) 

We note that the coordinate transformation given by (4.1) is a canonical transformation 
since the symplectic 2- form íl remains invariant: 

Q = dqi A dp\ + dq-i A dp 2 = dz A dp + dz A dp. (4-10) 



V. WIGNER FUNCTIONS OF HIGHER STATES 

Observing that the *-product (4.3) can be expressed in terms of a, d, b, b, explicit ex- 
pressions of the Wigner functions and other related computations can be carried out in an 
easier and more suggestive way To this end from Eqs. (4.4) and their complex conjugates 
we obtain 

z = -^J=(a + b), p = ^l{a-b), (5.1b) 
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and 



i y/2 

d z = -=(<9 a + <%), d p = (d à -d b ), (5.2a) 

0,= + a p = ^-(d a -d t ). (5.2b) 

7 a/2 mw7 



Then (4.3) and (4.7) can be written as 

* = exp[^(dada + Hfeis - dada - dhdb), (5.3a) 
H = -Tiu){aa + bb). (5.3b) 

By inserting (5.3b) and (4.8) into (3.22) and then by defining 

Wn = —a n * e~ 2aa * a n , (5.4a) 
n\ 

wi = ^b l *e- 2bi *b l , (5.4b) 

W n i can be cast to the following factorized form 

W n i = 4w n (a, a)wi(b, b). (5.5) 
By virtue of (5.3a) one can easily verify that 

n 1 1 

e - 2aS *a" = E(-ïï) fc Tïï(^ e " 2aS )( 9 a° n ) = 2 n a n e- 2aS , (5.6) 
fc=o 2 k - 

where (£) being a binomial number we have used the identity J2k=o(k) = 2 n - On substituting 
(5.6) into (5.4a) we obtain 

Oli 

Wn = —a n * (a n e- 2a ~ a ) 
n\ 

= —Aü - \d a ) n {a n e- 2aa ) = e- 2aà P n (2aa), (5.7) 

where we have defined 

P„(«) = ^e«(l-^)> B e- w ). (5.8) 
In the appendix we will prové the following operator identity 

(1 - ^-) n (u n e- u ) = (-l) n e u (^-) n (u n e- 2u ). (5.9) 

On multiplying both sides of these identity by e u / n\ and by recalling the Rodriguez formula 
of the Laguerre polynomials (see the appendix) 
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2u Jn 

Lni2u) = ^^r> ne ~ 2u ^ (5 - 10) 

we see that 

P n {u) = {-l) n L n {2u). (5.11) 

Therefore, from (5.7) and (5.11) the explícit functional form of diagonal Wigner functions 
(5.5) are obtained as follows 

W nl = {-l) n+l 4L n {4aa)L l {4bb)e^ aa+bi l (5.12) 



VI. SYMMETRY PROPERTIES 

Some transformation and symmetry properties of Wigner functions immediately follow 
from their factorized forms given by (5.12) and from following relations 

- HL 1 I M ■ -|2 

aa = — - = — — \p + iKz\ 
nuo ZTík 

[(pi + nq 2 ) 2 + (p 2 ~ tgi) 2 ], (6.1a) 



4Hk 

bb= — + — = 1 \p - i K z\ 2 

' [( P i-Kq 2 ) 2 + (p 2 + K qi ) 2 ], (6.1b) 



4Hk 

where k = muj/2. As is apparent from the first equalities of (6.1a) and (6.1b), the vàlues of 
W n i depend on the level sets of H L /huo and of J/U which are classical constants of motion 
and define some surfaces in R 4 . Rewriting (5.12) as 

W -„ 1 = (-ir4L„(4§)L I (4^ + 4§)e-^-&. (6.2) 

makes it possible to read out the vàlues of W n is on various level sets very easily. In particular, 
since L n (0) = 1, at the phase-space origin which corresponds to Hl = = J we have 

W^(0,0) = 4(-ir +í . (6.3) 

The following relations now easily follow from (5.12) and (6.1), 

WnK-q>p) = w ln {q, p) 

= W ní (q,-p), (6.4a) 
W nl (-q,-p) =W ní (q,p), (6.4b) 
W n i{q 2 ,qi,p 2 ,Pi) = W ln {q,p). (6.4c) 

(6.4a) shows that under the space inversion (q — > — q) and time reversal (p — > — p) transfor- 
mations all Wigner functions transform in the same way In particular, Wigner functions for 
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which n = l possess both space inversion and time reversal symmetries separately and they 
are positive valued. Eq. (6.4b) directly results from (6.4a) and it implies that all Wigner 
functions are invariant under the phase-space parity transformation (q, p) — > (— q, — p). 
The transformation property given by (6.4c) will be helpful in deriving marginal probability 
densities (see Sec. X). We should note that although (6.1a) and (6.1b) have some transla- 
tional invariance properties they do not lead to a translational symmetry of W n \. Namely, 
the following translational invariance of Hl. 

H L (q + c, p - kc*) = H L (q, p), 

where c = (ci,c 2 ),c* = (c 2 , — Ci), is not compatible with that of (6.1b). 

Due to the fact that W n i depend on the products aa, bb, a larger class of symmetries can 
immediately be identified. Indeed, it is evident that under transformations given by 



(a, ò, a, b) — > (ue^a, ve ir> b, u 



-ni 



b), 



(6.5) 



where £,77 and u 7^ 0,v 7^ are some real parameters, all W n \ remain invariant. Obviously, 
this transformation also leaves and J invariant. To express this transformation in terms 
of (q, p) coordinates, we define two column matrices x and y by x r = (a,b,a,b),y T = 
(gi, g2,Pi,í>2), where x T denotes the transpose of x. We then introduce two nonsingular 
matrices 

\ 



A 



(ue* 













u 



( 



B = ^ 



1 
1 

K 



1 

—K 






—i 
1 

K 

in 









-vq 



l 
1 

— K 

in ) 



i 1 



\ —IK —ÍK 

having the determinants detA = l,detB = 
transformation (6.5) and relation between y and x simply read as 

x' = Ax, y = £>x. 



(6.6a) 



(6.6b) 



= — h 2 . In this matrix notation the 



(6.7) 



In writing y = _Bx we made use of Eqs. (4.1) and (5.1). Since y' = 5x' = BAx. we deduce, 
from (6.7) that the transformation (6.5) can be written, in the canonical coordinates, as 



y 

c 

where y' T = (q[, <&,tfi,lQ and 



Cy, 

BAB- 1 



(6.8) 




M 




N 



s+ c_ \ 
-c 

c± = u + cos £ ± v + cos 77 + i (u- sin £ ± V- sin 77) , 
s± = u + sin^ ± v + sin 77 — i(u- cos^ ± t>_ cos 77), 

V 



(6.9) 
(6.10) 



u± 



u ± — , 

u 



v± 
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In view of detC = detA = 1 and 



NN T + MM T = 161, (6.11a) 
NM T — MN T = 0, (6.11b) 

where 1 and denote, respectively, 2x2 unit and zero matrices, it is straightforward to 
check that C is a symplectic matrix. That is, 

Jo= (-1 o)' 

being the Standard symplectic matrix, C satisfies the symplectic condition CJqC t = Jo- 
Hence C belongs to the symplectic group Sp(A,C). Evidently, W n i(Cïy) = W n i(y) = 
W n i(C2y) imply that W^CiC^V) = W n i(y). In other words, if C\ and C 2 matrices leave 
W n i(y) invariant so does CiC^ 1 . This means that the transformations (6.8) form a four 
parameter subgroup of Sp(4, C). 

Finally in this section we note that the volume form and symplectic 2-form in (a, b, a, b) 
coordinates take the form 

dV = dq\ A dqi A dp\ A dp2 = h 2 da A da A db A db, (6.12a) 
Q = dqi A dpi + dq 2 A dp 2 = ih(da A da + db A db). (6.12b) 

Then, the symplectic property of the transformation (6.5) can directly be observed from the 
second equality of (6.12b). 



VII. GENÈRIC PROPERTIES 



In this section we shall exhibit some genèric properties of diagonal Wigner functions, 
such as the reality, *-projection, normalization, and orthogonality properties, within au- 
tonomous framework of deformation quantization. Marginal probability densities and their 
bàsic properties will be taken up in section X. As is apparent from (2.3), (3.15) (or (3.22)) 
and (4.8) the reality condition W n \ = W n \ of diagonal Wigner functions is guaranteed by 
construction from the beginning. On the other hand, the following form of the so-called 
*-projection property 

W n i * W n 'i' OC 5 nn ,6 W W n i (7.1) 

can be easily inferred from (3.18) by taking the *-multiplication of them from the left and 
then from the right by W n >ii. We start by deriving the exact form of this relation for the 
Wigner function of Landau levels. 
Noting that, by (5.3a) 

k 2 

(aa) k+1 = (aa) k * aa + —(aa) k ~ l (7.2) 
one can verify (aa) k * W — (k\/2 k )W by induction and 
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-2taa 



* w ° = YTt Wo - 



(7.3) 



This is a special case of e~ 2taa * e" 2saa = (1 + ís) _1 exp[— 2(í + s)(l + ís) _1 aa] which can be 
proved by other means (see Eq. 34 in [19]). Hence, we obtain the *-projection property of 
Wq as follows 1 



'o * vvo 

By virtue of (3.13) we have 



Wn * Wq = Ae~ 2aa * e~ 2bb *W n = Wq. 



a n ★ a n ' = a"" 1 * a n '~ x * (N a + rí) 

= a n - 2 * a n '- 2 * (N a + rí - 1) * (N a + rí). 



This leads us to, by induction 



a n *a n 



a n '~ n -k Uf=n'-n+l( N a + j)*, for v! > 71, 

a n ~ n ' *Ilti(N a +j) in for rí<n, 
Il]=i(N a +j)*, for rí = n, 



where 



(7.4) 



(7.5) 



(7.6) 



U(N a + j)* = {N a + i)*{N a + i + l)*---*{N a + rí). 

j=i 

Similar relations hold for b l -kb 1 ' . On the other hand since 



(7.7) 



Wq -k N a -k Wq = = Wq -k N^-k Wq, 
Wq * é * Wq = = Wq * CL k * W , 



we see that 



I nn ifl> = W * a ra * a n ' * fe' * 6'' * W^o 

= n!Z!<J nn /<Jij/W * W = n\l\5 nn i8ü'Wo. 



Then, from (3.22) and (7.8) we obtain 



w nl *w n . 



1 



— a n *b l * I nn , :ll , k a n ' ★ 6 r = (W^'W^, 



n\l\rí\l'\ 



(7- 



(7.9) 



for all diagonal Wigner functions of Landau levels 

By virtue of (2.2), (3.22), (4.9) and (7.6) we have the following normalization property 



[24] a misleading factor 4/e 2 that appears at the end of Eq. (7.4) and in Eqs. (7.8), (7.9), 
(7.11) should be removed. 
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/ W nl dV = 4t7 / a n * b l * W ★ a n ★ ò'dV, 
Jr 4 n!í! Jr4 

= 4tt / a" * o n * * &' * WbdV, 
nu! ÍRi 

= / w dy = /i 2 . (7.10) 

JR 4 

Finally, the *-orthogonality directly follows from (2.2), (7.9) and (7.10): 
/ w nl * W n , v dV = í W nl W n/l ,dV, 

ÍR 4 iR 4 

= S nn .8 w í W nl dV = h 2 5 nn ,5 w . (7.11) 
Jr 4 

As is apparent, many properties of W n \ are guaranteed by the corresponding properties of 

Un. 

VIII. GENERATING FUNCTION FOR ALL WIGNER FUNCTIONS 

Let us define the phase-space functions 

G x [a,a) = e ai ~ a *e~ 2 ~ aa *e Pia , (8.1a) 
G 2 (b,b) = e Q2 ~ b * e- 2ib * (8.1b) 

where a. = (ai,a 2 ),(3 = (Pi, (3 2 ) are some parameters. Using the shorthands G\ = Gi(a,à) 
and G 2 = G 2 (b,b) we define and compute w niïl2 (a, à), wi^b, b) as follows 



w nin2 (a,a) = yJl/n^.n^d^ld^G^a^o^, 

= ^1/mWa 111 * e~ 2aa * a n \ (8.2a) 
w hh (b,b) = ^l/h\k\d l ^d%G 2 \ a2=G=P2 , 



'l/h\l 2 Ü h *e- 2bb *b h . (8.2b) 
We then define the generating function 

G = G 1 (a,a)G 2 (b, b), (8.3) 
from which all (diagonal and off-diagonal) Wigner functions can be generated as 

= 4iü nina (a,a)«; JlÍ2 (6,6), (8.4) 

In the case of n-y = n 2 = n, l\ — l 2 — l, (8.2) and (8.4) coincide with (5.4) and (5.5). 
Making use of 



e a í a * g-2aa _ e aia e ~ ± d a d a e -2aa _ g 2a(ai-a) ^ 
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one can easily verify that 

G 1 = e -°i/3l e 2(«iS+/3 1 a) e -2aa ; (g ^ 

and a similar relation for G 2 . Hence 

Q — e -a·f3 e 2(a 1 a+/3ia+a 2 b+/32b) e -^- ^ (8.7) 

where ol • (3 = a\f3\ + a 2 f3 2 and we have used (5.3b). 
Now let us consider the following expansion 

g-a^Ka+zM = g ^(2a) fe (l - |l) fc e 2/3ia . (8.8) 

fc=0 ' 

Making use of the generating function of the generalized Laguerre polynomials [28,29] 

00 

(i + y ) k e-*y = '£ L n- n (x)y n , (8-9) 

n=0 

which holds for \y\ < 1, we rewrite (8.8) as 

00 k 

e -ai0i e 2(aia+/ïio) = ^(2a)*- n (-/3i) n LÍ;- n (4aa). (8.10) 

k,n=0 

The condition |y| < 1 impose the constraint 2\a\ > \/3i\ for Eq. (8.10). Therefore, from (8.6) 
we have 

00 k 

Gi = e -2a~a ^ -±(2a) k - n {-^) n L k n - n {Aaa), (8.11a) 

k,n=0 fc - 
_ 00 k 

G 2 = e~ 2bb ]T T7 (26) fc - n (-/5 2 )"Lf t - ri (4òò), (8.11b) 

k,n=0 



and, by (8.2) 



'n 2 \ 



w nin2 (a,a) = W^-(-l) n2 (2a) ni - n2 L™r n2 (4aa)e-^, (8.12a) 



w hh (b,b) = J^(-iy-(2b) h ^L\l' l ^Abb)e- 2bl . (8.12b) 
V n- 

On substituting (8.12) into (8.4) we get 



W nin2hl2 = J^(-ir +l *(2ar- n >(2b) h - 1 *^ (8.13) 

We should note that when superscript of a Laguerre polynomials is a negative integer the 
following formula must be used 

L-\x) = (-x) k ^^L k _ k (x), (8.14) 

where k is positive integer. In such a case the resulting expression can also be obtained by 
starting with the following equivalent expansion of the left hand side of (8.8): 

00 R k ni 
e -ai/3i e 2(aia+/3ia) _ tLL[2a) k (l - — ) k e 2aiÜ . 

k=o k } - 2a 
In any case (8.13) coincides with (5.12) for n\ = n 2 = n and l x = l 2 = l. 
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IX. PHASE-SPACE COHERENT STATES AND GENERATING FUNCTIONS 
FOR MARGINAL PROBABILITY DENSITIES 



In this section we shall exhibit two important properties of the generating function G. 
The first is that G can be interpreted as a phase-space coherent state of Landau levels. The 
second property is that the integrated forms of G over some phase-space planes serve as 
generating functions for marginal probability densities (distributions) on these planes. 

For the first property we compute 

a * G = (a + ^d a )G = a±G, (9.1a) 
b*G= (b+-d- b )G = a 2 G, (9.1b) 

which imply that G is a left (complex) *-eigenfunction of the annihilation functions. On the 
other hand, since 

G*à = (a+ ^d a )G = fcG, (9.2a) 

G*b = (b+^d b )G = (3 2 G, (9.2b) 

G is at the same time the right *-eigenfunction of the creation functions. For these reasons 
G behaves as a left/right coherent state. Moreover 



a-kG = 


(a- 




(2a 


-(h)G = d ai G, 


(9.3a) 


b*G = 


(5- 


\d b )G = 


(26- 


-(3 2 )G = d a2 G, 


(9.3b) 


G-ka = 


(a- 




(2a 


- oíi)G = dfcG, 


(9.3c) 


G*b = 


(b~ 




(26- 


- &2)G = dp 2 G. 


(9.3d) 



provide us with the phase-space Bopp realization of left (right) *-actions of a, b (a, 6). To 
make G real it is sufficient to take a — (5. In such a case G corresponds to the Glauber- 
Perelomov Standard coherent state [25]. 

To exhibit the second property of G, we define the generating function for the marginal 
probability densities in the gig 2 -plane as follows 

/oo /*oo 
/ Gd Pl dp 2 . (9.4) 
-oo J — oo 

Then the marginal probability density P n i(qi,q2) in the gig2-plane can be generated as 

Pni(qi,q2) = ■^(d ai dp 1 ) n ( y d a2 dp 2 yM a p( y q 1 ,q 2 )\ a =o=f3- (9.5) 
Combining (9.4) and (9.5) and comparing the result with (8.4) we obtain 
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/OO Í-OO 
/ W n idp 1 dp 2 , (9.6) 
-oo J — oo 

provided that we can change the order of derivatives with respect to parameters with in- 
tegration on phase-space coordinates. The relation (9.6) justifies our naming M a p as the 
generating functions for marginal probability densities (see also (11.2)). 

Similarly, the generating function for the marginal probability density in the pijvplane 
can be defined as 

/oo roo 
\ Gd qi dq 2 , (9.7) 
-oo J — oo 

and the marginal probability density P n i(p 1 ,p 2 ) in the pip 2 -plane can be generated as 

Pm(pi,P2) = ^(d ai d^ n (d a2 d f32 ) l M a(3 (p 1 ,p 2 )\ a=0=f3 . (9.8) 

Generalizing these observations we can define generating functions M Q( g(gj,pfc); j, k — 1, 2 in 
order to generate the marginal densities P n i{Qj,Pk) on the g^-planes. 



X. MARGINAL PROBABILITY DENSITIES OF LANDAU LEVELS 

To ease the calculations of this section we shall use the dimensionless coordinates 

Z = -(qi+iq2), Z = -(q 1 -iq 2 ), (10.1a) 
7 7 

P= — ÍP1+ÍP2), P = —(P!-Íp 2 ), (10.1b) 

and define the dimensionless quantities 

p 2 = ZZ=\{ql + ql), (10.2a) 

7 



PP = ^ 2 (PÏ+PÍ), (10.2b) 



I q = («1 + fc)Z - (a 2 + p x )Z, (10.2c) 
I p = (a 2 - /3x)P - (ax - fo)P. (10.2d) 



We then rewrite (8.7) as 



Substituting (10.3) into (9.4) we obtain 

M aP {q u q 2 ) = e-^e~^K x K 2 = 7r(-) 2 e^ Q a Q p , (10.4) 

7 



where 
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J — oo 



TrVaígèíai-aa+^i-A) 2 

7 



= r e -^ [ ^ + ^ (ai+a2 - ft - /32)p2] ^ 2 

J — oo 

= TrVa^g-èíai+aa-ft-A) 2 

7 

are obtained by the Standard complex integration methods and 

Q a = e -«l«2+'(«l^«22) ] = e -W2~l{^Z-P 2 Z) _ ^ 

When Eq. (10.4) is substituted into (9.5), in terms of 

^ní = ^i^ia^aUi^aa, ^ní = % #g 2 <2/3 1 =0=/3 2 > (10.6) 



we obtain 



47T 7) q 

^(gi,g2)^(-)V"Ji4. (io.7) 

i t.i . y 



Similar to (8.8) and (8.10), we can expand (10.5) as follows 



= £ Ç^'- fc a{a 2 fc Lr fe (p 2 ), (10.8a) 

j=o -? ! z 

= E t^^-^^Lr fc (p 2 ). (10.8b) 

j,fc=0 

Then, by using these in (10.6) we have 

Jl = i n+l l\Z n - l Lr l (p 2 ), J 2 nl = (-t) n+l l\Z n - l Lr l (p 2 ), (10.9) 

and, by (10.7) 

Pm(qi,q2) = An l ^fp^-%-f[Lr\p 2 )?. (10.10) 

The generating functions M a p(pi,p 2 ), M a p(qi,pi), M a p(q 2 ,P2) and the corresponding 

marginal probability densities can be calculated in a similar way These are all calcu- 
lated and presented together in the Tables I and II. But the calculations for the qip 2 and 
<?2Pi-plaries are to be presented as they are a bit different. 

From (10.3) we have 
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/OO Í-OO 
/ Gdq 2 dp 1 , 
-oo J — oo 

= 7r ^ e i(r2+T2) e I[( Ql +ir_) 2 + (/3i-ir_) 2 + (a 2 -ir + ) 2 + (/32+ír + ) 2 ] ) (10.11) 

where t± = (mugí ± 2p 2 ) /mu^. We now consider the identity 



= {-iTd n T _e-^- = (-L)» e -^lf B (^). (10.12) 

In passing to the last equality we have used the definition of the Hermite polynomials 

H n (x) = (-l) n e x2 d^e- x \ (10.13) 



Then by using (10.12), (10.13) and the other three relations obtained for derivatives with 
respect to a 2 , Pi and f3 2 hi 

Pni(qi,P2) = ^(d a M n (d a M l M aP ( qi , P2 )\ a=0=l3 , (10.14) 

we arrive at 

As a common property, all marginal probability densities of Landau levels are positive 
on the corresponding phase-space plane. They determine the wavefunctions on these planes, 
up to some unitary phase factors. Axial symmetry on the respective planes of the first 
four marginal probability densities are obvious from their explicit expressions given in Table 
II. Also in functional dependence P n i(qi,P2) and P n i(q2,Pi) differ from others. Marginal 
probability densities along qj, or pj, or along some oblique phase-space directions can be 
computed as well. Finally we note that P n i(<l2,Pi) presented in the last row of the Table II 
easily follows from (10.15) by virtue of (6.4c). 



XI. REMARKS 

In this section some remarks that remain outside of our main goals are collected together. 



A. Comparison with Definitions 

For justification we would like to compare (4.8) and (10.10) with expressions obtainable 
from the definition of the Wigner functions. In the symmetric gauge (3.20) and in terms of 
plane polar coordinates (r, 0) the normalized wavefunctions and the corresponding energy 
eigenvalues of the Landau levels are given by 

ú .( r fí) = -. I Url \j\ e i36 - P 2 /2 L \j\ / 2X 

1 I7I-7) 

E nr ,=hu(n r + -+ l -^ Y ^), (11.1) 
n r = 0,1,2,..., j = 0,±l,±2,... 
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where r = 7p = (qf + q^) 1 ^ 2 and n r ,j stand for the radial and angular momentum quantum 
numbers. Since in (10.10) it is supposed that n > l, we can make the identifications \j\ = n—l 
and n r + \j\ = n, which imply n r = l. We then have, from (10.10) and (11.1) 

Pm(quq2) = h 2 \^ nrJ (r,e)\ 2 . (11.2) 

This expected result confirms Eq. (10.10) and our definition of marginal probability densities 
in section IX. On the other hand, by using the normalized ground state wavefunction V'o(q) — 
^ooÍT-,^) — 77T -1 ^ 2 exp(— p 2 /2) in the definition 



r l — l 

W ° = 7r2 ^°^ q + 2 y ^°^ q ~ 2 y ^ e ^ dyidy2 ' ^ 1L3 ^ 

we get 



W = ^r e-*<»«*»r dyi r e -^^ dy2 . (11 .4) 

Noting that, each integral in (11.4) is equal to 2^-n 1 / 2 we obtain W = Aexp(—AH /huj) 
which is the same as (4.8). However, as we have mentioned in the introduction, in attempt 
to find other confirmat ions the resulting integrals do not seem to be so easy to cope with. 

B. Unitary Transformations in a Phase-Space 

Let us consider a unitary similarity transformation generated by the phase-space function 
U = C/(q, p) and its *-inverse U~ x = U: 

U + U- 1 = U~ 1 *U= 1. (11.5) 

When such a transformation applied to a *-eigenvalue equation H *W\ = W\ * H = E\W\, 
the transformed functions 

H' = U * H *U~\ W^U-kWx-kU- 1 , (11.6) 

obey the same *-eigenvalue equation 

H'*W' X = W' X *H' = E X W' X . (11.7) 

In view of (2.3), (7.9) and (11.5) it is easy to verify that 

W' x = U -kW\ -k U -1 , (11.8a) 
W' x * W' x , - kW' x 5xx' = U*(W X * W x > - kW x 5xw) * U~\ (11.8b) 

where k is a constant. That is, the reality and *-projection properties of W\ are preserved 
under the unitary similarity transformations. Moreover, as is evident from 

/ WídV = í U * W x * U~ 1 dV 

JR2D A J R 2D 

= í W x *(U- 1 *U)dV= í W x dV, (11.9) 

JB?D JR2D 

normalization properties of Wx is also preserved. In (11.9) we have used (2.2) and (11.5). 
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C. Gauge Transformat ions in a Phase-Space 



Now let us consider 

U q = e^ x , U' 1 = V q = e-^ x , (11.10) 

where U q = U(q) and x = x(q) is a rea l valued function of the coordinates. Obviously 
U q -k U~ x = U q U q l = 1. By using (1.3) with D = 2 we also have 

tf ? *p*£C 1 =P--v ï x, (ii.ii) 



c 



which implies, for the Landau Hamiltonian (3.1) 



H' L = ^P^(P- ~ c A ) * U, 1 ] ■ P 9 * (P " ?A) * O 

1 r__ 9/» r- ,i2 



= ^[P-^ A + V * < 1L12 > 

where H' L = U q -k * C/" 1 and V g represents two dimensional gradient operator in çi,Ç2 
variables. Eqs. (11.11) and (11.12) show that the unitary transformation by C/ g amounts 
to the so called gauge transformations of the second kind in p and Hl- Hence, without 
changing their genèric properties the Wigner functions of Landau levels in any gauge x can 
be computed directly from 

W' nin2hh =U q * W nin2hh * U q \ (11.13) 

where W / ni „ 2 ; 1 ; 2 are given by (8.13). In that case the transformation in (11.13) may be 
considered as a first kind gauge transformation. We should emphasize that the two forms 
of *-product given by (1.3) and (5.3a) can equally well be used in (11.13) provided that the 
arguments of functions are properly written. 



D. Implications for Wigner Functions 

By using (1.3) we compute the following useful relation in an arbitrary dimension 
/i(q)*e-*™*/ 2 (q) = [e-^ p+ ^) /l(q) ] * / 2 ( q ), 

= /i(q+^y)[e-* y " p */ 2 (q)], 

= /i(q + ^y)/ 2 (q-iy)e-^ p , (n.u) 

where V g denotes D-dimensional gradient operator in Çj variables. Taking f\ = U q and 
f 2 = Ug 1 in (11.14) yields 

C/ ï *e - * y,p + = e i ^ x(q+ 5y) e - í ^ x(q -5y) e -s y · p (11.15) 
This implies, for W' XlX2 = U q -k W Xl \ 2 *U q l , that 
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= /^ifq + ^fq-^^^n!/), (n.16) 

where W\ 1 \ 2 is a off-diagonal Wigner function defined by (1.1) and ^.(q) = e^ x ^xjj\ j (q) = 
Uqipx^q). If we change our notation as W^ x ^ = W\ 1 \ 2 , Eq. (11.16) implies the following 
important relation for gauge transformation of the Wigner function 

= U + W^+U- 1 = W^, = W Uq ^, Uq ^. (11.17) 

Finally we should note that, for f\ = ipi and f 2 = ïp 2 Eq. (11.14) suggests the following 
redefinition of Wigner function 

W^ 2 = í M^*e-r^^ 2 ( q )dV(y), (11.18) 

which enables us to intèrpret the map W as a "sesquilinear *-Fourier transformation" . 

XII. CONCLUSION 

In this paper we have presented a problem that can, with all of its features, thoroughly 
be considered in the framework of deformation quantization without making use of other 
quantization schemes. As a summary, we have obtained Wigner functions of Landau levels by 
solving a two-sided *-eigenvalue equation, specified a large class of its transformation and 
symmetry properties as well as established its genèric properties in the same framework. 
Off-diagonal Wigner functions and marginal probability densities are generated by means of 
generat ing functions. 

With its distinguishable and intriguing properties, Wigner function provides the most 
complete description and offers convenient interpretation of the result of several recent ex- 
periments in atòmic òptics, molecular physics and signal processing [7,8,26]. A bàsic aim 
of all these experiments is to probe the fundamental structure and predictions of quantum 
mechanics in a new way, by observing non-classical behaviors of Wigner function [27]. For 
this purpose marginal probability densities play equally important role since Wigner func- 
tion are reconstructed from measured densities along various phase-space directions . Now 
explicit expressions of Wigner functions and all marginal probability densities of Landau 
levels at hand we do expect that the results of this paper will shed more light on future 
experiments in which Landau levels are realized and on their interpretations. 
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APPENDIX A: THE OPERATOR IDENTITY (5.9) 

In this appendix we will prové the following operator identity (see Eq.(5.9)) by induction: 

(1 - T) n (x n e' x ) = (-l) n e x T n (x n e- 2x ), (Al) 

where T = d/dx and n is a positive integer. Then, we will present a finite sum formula for 
the generalized Laguerre polynomials that follows from (Al). 

Eq. (Al) holds for n = 1 (and obviously for n — 0). Now we suppose (Al) for n — 1 and 
try to prové it for n. Defining — (1 — T) k (x k e~ x ), the inductive hypothesis is 

I n _ x = (-l) n ' 1 e x T n - 1 {x n - 1 e- 2x ). (A2) 

Then we have 

I n = (1 - T) n [x(x n ~ 1 e~ x )] = [x(l — T) — n]/ n _!, 

= (-l)V[ir + nT n - 1 ](x n - 1 e- 2x ), (A3) 
= (-l) n e ie 7 ,n (a; n e- 2iB ). 

In passing to the second equality we made use of [(1 — T) n ,x] = —n(l — T) n_1 , where [, ] 
stands for the usual commutator. The second line of (A3) easily follows from (A2) and in 
passing to the last line we have used 

{xT n + nT"' 1 )/ = T n {xf), (A4) 

where / is an arbitrary differentiable function of x. Eq. (A4) easily results from the 
commutation relation [T fc ,x] = kT k ~ l , or equivalently from the Leibniz rule T n (xf) = 
Ek=o(k)( Tk x)T n ~ k f. We should note that (Al) can be generalized as 

(1 - T) n {fe x ) = {-l) n e x T n f, (A5) 

which can be proved more easily, again, by induction on n. (Al) corresponds to choice 
/ = x n e~ 2x in (A5). 

It is worth mentioning that, in view of the binomial expansion 

(1 - T) n = è(7)(-l) n " J 'T n - J ', (A6) 

3=0 

and by recalling the Rodrigues formula for generalized Laguerre polynomials ( [28], p. 241) 

L a n (x) = —x- a T n (x n+a e- x ), (A7) 

71/. 

from (Al) and (A7) we obtain the following finite sum formula 

Ln(2x) = £ L-Ll^x) = X /./(•<•!/.;, ; í.r). (A8) 

3=0 3' 3=0 

where LJ J (x) = (—xY/j\. Since L n = L°, from (A7) we get (5.10) for a = and x = 2u. 
We could not find (A8) in the related classical references [28,29]. 
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TABLES 



TABLE I. Generating functions for the marginal probability densities for Wigner functions in 
various phase-space planes. 



M af s(p 1 ,p 2 ) = 7 r7 2 e - 4 Ç 2 e «i«2+2(QiP-a 2 P) e /3i/32+2(/3iP-/32P) ) p = P1+ÍP2 



M a p( qi , P 2) = vr ^ e |(-++--) e è[(°i+--) 2 +(/ 3 i---) 2 +(«2— +) 2 +(^+-+) 2 ], r± = ^ ± 
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TABLE II. Marginal probability densities for Wigner functions in various phase-space planes. 
In this table we use the abbreviations N n i = Airl\/n\ and N' nl = Air / n\l\2 n+l . 



Prabuto) = N nl ^?p^- l )e-?[Lr\p 2 )] 2 , 


„2 _ ll+lï 
P ^ 


Pnl(pi,P2) = iV„a 2 (2C) 2(n - /) e- 4ç2 [L[ l -'(4C 2 )] 2 , 


z-2 _ ■y'Hp'í+PÍ) 
* au 2 


Pm(qi,Pi) = N nl h^ n - l) e-iï[Lr l (pÍ)? , 


M? = £ + £ri 


Pnl(q2,P2) = N nl h$ n - l) e-Ú[L?- l (j4)? , 




Pm(qi,P2) = Kfie-m+^H^H^f . 


T-, _ 91 _|_ 7P2 

T ± — 7 ^ h 


PnM2.P1) = Kfie-^+^HX^H^)? , 


t-/ _ 92 1 7P1 
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